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We explain the PAMELA positron excess and the PPB-BETS/ATIC e+ + e" data using a simple 
two component dark matter model (2DM). The two particle species in the dark matter sector are 
assumed to be in thermal equilibrium in the early universe. While one particle is stable and is the 
present day dark matter, the second one is metastable and decays after the universe is 10~* s old. 
In this model it is simple to accommodate the large boost factors required to explain the PAMELA 
positron excess without the need for large spikes in the local dark matter density. We provide the 
constraints on the parameters of the model and comment on possible signals at future colliders. 



INTRODUCTION 



Other than its gravitational interactions, the proper- 
ties of dark matter are largely unknown - only lower limits 
on the mass and upper limits on the couplings have so far 
been obtained. It has been suggested that the prelimi- 
nary results from The Payload for Antimatter Matter 
Exploration and Light-nuclei Astrophysics (PAMELA) 
experiment may represent a breakthrough in this situ- 
ation PAMELA sees a larger positron fraction in 
the cosmic ray flux at 10 - 80 GeV than one expects 
in the galactic environment, in agreement with previous 
hints from the HEAT and AMS-01 experiments [3, 
A possible explanation for this excess is that it comes 
from dark matter (DM) annihila ting in our own galaxy 
HHQIlJor from DM decay- 



ing with a d ecay time much loriger than the age of the 
universe IigI. [itL [3, [3, 113, HE 111] (or the combination 
of the two [22] )■ Astrophysical sources, such as a nearby 
pulsar, have also been su ggest ed as an explanation of the 
excess 



H, Si m m miii, m, H^. in this paper we 
will focus on the annihilating DM interpretation of the 
excess. 

The PAMELA experiment also detects an anti-proton 
cosmic ray spectrum which is compatible with the ex- 
pected galactic background, giving a constraint on dark 
matter annihilation modes. The authors of Ref. [5| find 
two explanations for the PAMELA data: i) a heavy DM 
particle of mass above about 10 TeV that annihilates 
predominantly to a W^W~ pair or ii) a DM particle 
which annihilates predominantly into SM leptons with 
no strong constraint on the DM mass. 

There are also hints of a positron excess in the energy 
range of 500 — 800 GeV from a long duration Polar Patrol 
Balloon (PPB-BETS)flight and from the ATIC-2 bal- 
loon experiment [sil . |33 ] , in agreement with the change 

33| . If one 



which annihilates into leptons [1, H, 0, [1] . 

One class of DM candidates are particles that have 
a weak-scale self-annihilation cross section at freeze-out 



{aAv)Y 3 X lQ-^^cw?/s. 



(1) 



The resulting thermal relic abundance is then close to 
what is required to explain dark matter. Such a can- 
didate does not, however, fit straightforwardly with the 
PAMELA observations. If the PAMELA excess flux of 
positrons, is due to DM annihilating in the milky 

way, then 



2 

DM' 



(2) 



where pdm is the local DM density, ttidm the DM particle 
mass and the average number of positrons produced 
in a single XX ~^ annihilation (usually Ne+ < 1 and 
at most a factor of few in more exotic models). The 
proportionality in ^ includes among others the effect 
of positron propagation in galactic medium. For masses 
above 100 GeV the PAMELA positron excess suggests a 
(aAv) which is larger than the value required at freeze- 
out {aAv)F [IjB, 7, 8]. The mismatch is parametrized by 
a parameter called the boost factor 



B = N, 



(3) 



where poM — 0.35 GeV/cm"^ is the average expected local 
DM density. We can have B > 1 if either {(jav) > {<jav) f 
or if there is a local DM over-density due to substruc- 
ture in the galactic distribution^. Analysis of the N- 
body dark matter only Via Lactea II simulation suggests 
there is a 1% probability of density fluctuations enhanc- 
ing the local annihilation rate by a factor S ^ 10, while 



of power-law seen by HESS above ~ 1 TeV 
tries to explain simultaneously both these data and the 
PAMELA data using DM annihilations, it seems neces- 
sary to consider a DM particle with mass around 1 TeV it as referring to the combined effect. 



^ Conventionally, "boost factor" refers only to this astrophysical 
origin of the cosmic ray enhancement, while we prefer to define 
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much smaller fluctuations were observed in the Aquar- 
ius Project simulations [3^|3l|. While substructure in 
simulations is not well understood, larger boost factors 
seem to require a different explanation - an enhanced 
DM annihilation cross section. The values of B which 
explain the positron data are roughly B ^ (10, 10"^, lO'^) 
for moM (100 GeV, 1 TeV, 10 TeV) respectively, with 
a variation of a factor of few depending on the predomi- 
nant annihilation channel 0, 0] . 

If the annihilation proceeds through the s— channel, 
then a A o: 1/v and {aAv)F at freeze-out is ~ {oav) to- 
day, giving no enhancement in ([3]). If the annihilation 
proceeds through the p— channel then ct>i (x w meaning 
that {(Jav) <C {aAv)F- In both cases the DM interpreta- 
tion of the PAMELA data for ttidm > 100 GeV seems to 
exclude a simple single component thermal relic explana- 
tion of the results. This is true also, if DM is composed of 
many stable components as in 37| , in which case B <\. 
In the case of co-annihilation, the boost factor is bounded 
by i? < iV|, where Nf is the number of co-annihilating 
flavors. Again, large boost factors are excluded unless 
Nf is very large. 

If, however, the s— channel annihilation is enhanced by 
Sommerfeld corrections, then [gav) ^ (o'At')F is possi- 
ble. For DM particles which interact with VF, Z bosons, 
these Sommerfeld corrections are present if the DM par- 
ticles are heavy, itidm ^ Airmw 1 9^ ~ 2 TeV. The annihi- 
lation cross section is then enhanced by jv and can be- 
come large for w ^ [s^ . Alternatively, the Sommerfeld 
enhancement can come from a new force in the DM sector 
where a force carrier possesses GeV mass 0, with very 
large Sommerfeld enhancements possible, if there exists 
a bound state very close to threshold [s^, [i^ . Enhance- 
ments are also possible if the annihilati on g oes through 
a resonance with mass of order 
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All these explanations require light new degrees of free- 
dom with mass below moM- We will address the other 
possibility, namely that the apparent enhancement comes 
from states that are heavier than tudm- As a working 
tool we will present a simple 2-component dark matter 
(2DM) model that can explain the PAMELA and balloon 
experiment data ^. 



TWO COMPONENT DARK MATTER 

In the 2DM model the dark sector is composed of two 
particle flavours. These two DM particles can be scalars, 
fermions or vectors. The first DM particle xi has mass 
mi and is stable. We assume that it is this particle which 
is the cold dark matter (GDM) relic still present today. 



^ The name is somewhat of a misnomer since only one of the two 
components is stable on cosmological timescales. 



responsible for galactic rotation curves and also, through 
self-annihilation, the PAMELA and balloon results. The 
second DM particle X2 has mass m2 > mi, is unstable 
and decays to xi- We also assume that both X2 and xi 
are in thermal equilibrium prior to freeze out. We will 
discuss under what conditions such a setup is able to ex- 
plain PAMELA and balloon data. Many decay modes of 
X2 are possible including a two body decay X2 2xi or 
multibody decays X2 ^ Xi + AT, where X are SM parti- 
cles. To retain generality we define iVdec as the average 
number of xi particles produced in a single X2 decay. 

The number density of each of the dark matter flavours 
evolves with time as follows 



dni 

~dt 

dn2 

~dt 



(4) 



where nicq(t) is the equilibrium number density given by 
the thermal Boltzmann distribution and r2 is the decay 
width for X2 NdccXi + ^ decay. It is useful to de- 
fine z = TTii/T and normalise number density to entropy 



density, Yi{z) 
tions become 

z dVi 



ni{z)/s{z 



Then the evolution equa- 
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2oq 



[crAiVi), while 



where the annihilation rates are Tai 

Yicq{z) = nioq{z)/s{z). 

An interesting limit to consider is F2 <^ ^Ai,2- The 
abundances then first settle into their thermal relic val- 
ues, Yi Y^^ '"^^-, after which X2 decays to xi (see Fig. 
[ij. In this limit we then have 



Yi(oo) = Yi- 



Th.rcl. 



+ NdecY2 



Th.rcl. 



(6) 



The number density of the dark matter relic is enhanced 
by the number density of metastable dark matter compo- 
nents. Using an approximate analytic solution for s-wave 
annihilation, where y^Th.rei. ^ l/(^crAiVi), we find 



Yi{^) 

y^Th.rcl. 



l + A^dcciifl-^^logi?), (7) 

^ Zi ' 



where ~ 20 is the freeze out value of z for the stable 
DM component, while 



R 



m-i {aAiVi) 

TO2 {crA2V2) ' 



(8) 



We can now see why it is possible to explain the 
PAMELA and balloon data using the 2DM model. If i? > 
1 we will have Yi (00) > y^Th.roi. ^ rpj^^ positron fiux excess 
measured by PAMELA and given by Yi(oo)^(cr^iUi) will 
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FIG. 1: The solution for Y^{z) (black solid line) and Y-i{z) 
(blue solid line) for mi = 1 TeV, m2 = 3 TeV, r2 = Vd~'^^ 
GeV, A'doc = 1 and three different values of R as denoted. 
Dashed lines denote the thermal relic values of Yi(_z). On 
the upper figure {(JA1V2) is held fixed to ((tau)cdm = 3 x 
10~^®cm^/s, so that without decay this would give correct 
DM relic density with X2 the DM particle. Through decay 
this is transfered to xi- For illustration we also show the 
lower figure, where (ctaiWi) is held fixed to 3 x 10~^^cm'^/s. 

therefore be larger than one would expect in the case of 
one-component DM model, where the positron excess is 
proportional to (^yTh.rci.~j {oax'^x]- This is shown on Fig. 
[TJ The larger the ratio the larger the enhancement of 
Yi(oo). Note also that as the value of r2 approaches the 
values of V a\,A2 there is a washout effect, as seen for the 
R = 10^ curve in Fig. [TJ If X2 decays too quickly, then 
the resulting xi particles may still have chance to anni- 
hilate with each other before xi freezes out completely. 
If we are not therefore in the limit where r2 <C ^ai,A2, 
the enhancement effect is lost. 

We also show in Fig. ^ the predicted boost factors 
in the limit of small r2 



B = 



((Tav)f 
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yTh. 
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(9) 



The first two ratios on the right-hand-side are 0(1), while 
the last ratio is given in ([7]). The quantities labeled CDM 
are the parameters for the usual WIMP cold dark matter 
scenario and 2DM corresponds to the parameters of our 
model. 

We next discuss the limits on the parameters of our 
model. In order to explain the PAMELA data one needs 
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FIG. 2: The boost factor _B as a function of mi and R, in the 
case where r2 is small enough that wash-out can be neglected. 
The framed numbers (n) labeling contours give boost factors 
as B = 10". 

mi to be larger than 50 — 100 GeV, with mi ~ 1 TeV 
suggested by the ATIC excess. For X2 to be able to decay 
into xi we need 7712 > mi. For a particular assumed 
annihilation channel XiXi ATsm, the PAMELA and 
balloon experiment data fix {(Tai'Vi) — B{mi){(7Av) f as 
a function of mi Q. The value of m2(cr^2'^2) then follows 
from Fig. [51 For large boost factors we have {(TA2V2) ^ 

{cTAlVl). 

To have significant boost factors, the washout effect 
should be small, giving an upper limit on r2 . This means 
that at freeze-out, zf", the second term in ([5]) is smaller 
than the first one (we approximate Y2{z) at late times 
with its freeze-out value) 



r2 < 



Tai{z[)Yi,^{z[) Tai{z[) 



^2cq(zf ) 



B 



(10) 



and the annihilation rate at freeze-out is approximately 

2 

r^i(zf)^3i/(zf)~0.1^, (11) 



mpi 



where in the last equality we used z[ ~ 20 and ~ 10^. 
Equations (jlO[) and (|lip imply an upper bound on r2 



r2 < lO^^GeV 



(!)■( 



mi 



ITeV/ 



V 



or 



7 , / B \ /lTeV\2 



(12) 



(13) 



Because of this, if X2 is produced in colliders at relativis- 
tic velocities, it will travel at least a few meters before 
decaying. 

There also exists a lower bound on r2 from nucleosyn- 
thesis. One can use the detailed results presented in [44 1 
to argue that so long as the lifetime of X2 ^ Xi + SM is 
less than around a second there will be no change to 
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light element abundances. If on the other hand only 
X2 XiXi is allowed, this usually does not affect abun- 
dances. The X2 decay time can then be very long, of 
order the cosmological time scale. The simplest scenario 
then is that by the time of structure formation the xi 
particles from the decay are non-relativistic so that they 
are cold dark matter. 

More interesting is the situation where the particles 
decay after the start of structure formation. In the 1980s 
a significant amount of thought went into the idea of de- 
caying dark matter and its effect on structure formation. 
The motivation at that time was the idea that the missing 
energy in the universe now usually subscribed to being 
dark energy was actually the relativistic decay products 
of dark matter 4^ 3, 43] . Such particles would not clus- 
ter below 100 Mpc and therefore would not contribute to 
the value of JIm ^ 1 measured on scales smaller than 100 
Mpc. This scenario is disfavoured by more recent obser- 
vations of II(z) and structure formation and would not 
anyway solve the problem of the excess positrons we are 
addressing in this work. However, work done at the time 
and more recently has shown that the decay of dark mat- 
ter into another dark matter species (cold or hot) and its 
heating in that decay has interesting effects on structure, 
puffing out dark halos making them more diffuse (isl . [4^ . 
Because of this some workers have suggested that decay- 
ing dark matter could have a bearing upon two apparent 
possible problems with the KCDM scenario, namely the 
cuspyhalo problem and the small scale power problem 
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l52l |. While these problems are controversial in 
that not everybody actually agrees if they exist, it is cer- 
tainly true that the possibility of dark matter decaying 
into another dark species has a rich phenomenology when 
it comes to structure formation. 



PARTICLE PHYSICS CONTEXT 

To recapitulate, we have found that the 2DM model 
with two DM components, xi that is stable, and X2 
that is metastable, can explain the enhanced annihila- 
tion cross section observed by PAMELA. The large boost 
factors observed are explained by the hierarchy 



mi {aAiVi) ^ mi A| ^ ^ 
m2 {aA2V2) m2 5I A^ 



(14) 



where we have denoted schematically the dependence of 
cross sections in the non-relativistic limit on couplings 
gi^2 and masses Ai,2 of exchanged particles. The hier- 
archy of annihilation cross sections {aAiVi) ^ 
can be obtained for instance, if the typical coupling gi 
in the first DM sector is larger than 32 of the second 
one (here the hierarchy need not to be very large, for 
instance even for a boost factor of 10^, 32 ^ O.lgi suf- 
fices). The other possibility is that the annihilation of xi 
proceeds through a heavier state than the annihilation of 
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Xi 


X2 


SM 


Z2 


-1 
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Z'^ 
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-1 


1 



TABLE I: The charges under ^2X^3. 



X2- For large boost factors a relatively large hierarchy 
is needed, though (for instance for a boost factor of 10^, 
Ai ^ 10~^A2). Another elegant possibility is that the 
annihilation of xi proceeds through and s-wave process, 
while the annihilation of X2 is p-wave suppressed. The 
annihilation cross section for X2 is then suppressed 
(see e.g. [i^), which for v ~ 0.05 at freeze out can 
lead to a boost factor of ~ lO'^ without any fine-tunings. 
This is easily realized in a concrete model, if xi is a 
Dirac fermion, while X2 is a Majorana fermion. Yet an- 
other possibility is, if m2 annihilation is phase space sup- 
pressed. For instance, if m2 and m\ are almost mass de- 
generate and X2X2 XiXi is the dominant annihilation 
channel for X2i while xi can also annihilate to SM. 

More severe is the hierarchy between r2, Eq. (fT^ . 
and the decay width F ^ g^mj ~ few GeV typical 
for a weakly coupled theory. One explanation would be 
that X2 carries an approximately conserved charge that 
suppresses its decay. A simple possibility is that Xi is 
charged under Z2 and X2 under a different ^2' while the 
SM is neutral under Z2 x Z'^. Then the annihilations 
XiXi XsM, X2X2 X are allowed, while xi,2 ^ 
X decays are not. If Z2 is broken at some high scale 
A ^ 1 TeV, X2 is metastable in agreement with the 2DM 
explanation of the PAMELA data. Another possibility 
which we do not pursue here is that X2 could be charged 
under a gauge group broken at a high scale. 

We next discuss the possible interactions xi and X2 
can have. For simplicity we focus on the case, where 
both xi and X2 are scalars. If xi,2 are singlets under 
SU{2)l X U(1)y, then the most general renormalizable 
Lagrangian invariant under Z2 x Z2 is 



= Ain + ClX?X2 - 

+ C3{H^H)xl 



(H^H)xl 



C2 



4 

C4X1 



C5X2- 



(15) 



The Xi and X2 thermalise through interactions with 
the SM higgs and through four scalar interactions 01^4,5 
in the dark sector. In the milky way xi annihilate 
into hh, WW, ZZ leading to leptonic and hadronic final 
states. Thus xi mass has to be large enough, mi ~ 10 
TeV so that xiXi ~^ W^W~ can explain both the 
e+/ (e++e~) excess as well as the absence oip/p signal by 
PAMELA (while the balloon experiments, if confirmed, 
would exclude this simple scenario) Q. It is possible to 
avoid this constraint by enlarging the dark sector, for in- 
stance by charging xi under an extra U(l) under which 
also the SM leptons are charged but not the quarks, while 
keeping C2 small enough. 
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To get large boost factors in this simple scenario the 
hierarchy |c2| » \ci,3 \ is needed, for instance, for B ~ 10'' 
a hierarchy \c2\ ^ lO^^jci^sl works. Also, in order to 
prevent "fast" X2 decay, the Z2 should not be broken by 
a X2 vev (i.e. (X2) = 0, which for instance is trivially 
true if all q > 0). These qualitative conclusions do not 
change even if xi and/or X2 are neutral components of 
some higher representation of SU{2)l x U{1)y- Even 
then the terms invariant under Z2 x Z2 will contain two 
Xi or two X2 fields. Renormalizable interactions can thus 
couple xi,2 only to the higgs or to W, Z. 

The dimension 5 operators relevant to the X2 decay 
that break Z2 but not Z2 are 

^(V;{1,75W)X2, '^{H^H)x2, \{H^Hfx2, 

(16) 

where ip are the SM fermions. The dimensionful parame- 
ter TO is most conveniently chosen to be ~ TO2 (the choice 
of TO only rescales the definition of A). Decays into all 
massive standard model particles are possible with par- 
tial decay widths 

T2 ^ ml/{l&nK^). (17) 

The last and the first operator in give decay widths 
that are additionally (w/to2)^ suppressed. From Eq. (fT2|) 
and (jl7[) we then have (taking toi ^ m2 = toi_2 for 
simplicity) 

1, / B \ i/2/ mi2 \i/2 , ^ 

A»10-GeV.(^) (^) . (18) 

This general limit applies also, if Xi,2 are fermions or 
vectors. 

If we assume that X2 decays with T2 < Is (so that BBN 
constraints are fulfilled regardless of the decay mode), 
this then implies an upper limit on A in (|17p 

A<5a0-G.v(^)"', (19) 

Note that the allowed range includes the see- 

saw scale < 10^^ GeV, an intriguing possibility in view 
of the leptonic-only signal in PAMELA. 

We also briefly comment on the case where xi and/or 
X2 are fermions, while leaving detailed analysis for fu- 
ture. In order to thermalise with the SM, the fermionic 
Xi,2 need to be either charged under the SM gauge groups 
or couple to a hidden sector that then mediates with the 
SM. The simplest case is that xi and/or X2 are weak dou- 
blets so that the neutral components are massive "dark 
neutrinos" (some dark multiplet must also exist to can- 
cel anomalies). Then xi and/or X2 have masses in the 
10 TeV range as in the scalar case. Again, this can be 
avoided by enlarging the dark matter sector. 



What would be the signatures of a 2DM model in a col- 
lider experiment? If X2 is produced near threshold, it can 
decay in the detector and may be observed directly. Let 
us be more specific with a few illustrative examples. If 
Xi is a singlet scalar and X2 a weak doublet fermion, then 
the decays X2 Xi^^ are possible. If on the other 

hand xi is a weak doublet fermion, while X2 is a singlet 
scalar, the decays X2 ^ 2xi, X2 ^ XiXi ^ -'^SM and 
X2 Xi'^j Xi^^ I 7 99 are possible. If xi,2 are both 
weak doublet fermions, then the possible two body decays 
are X2 ^ Xi^i^Xi^: X2 ^ vh.vZ, X2 ^ xtW~,l+W~. 
How challenging the experimental search for X2 may be 
will depend on the actual masses and branching ratios. 



DISCUSSION AND CONCLUSIONS 

The proposed 2DM explanation of PAMELA/ATIC 
anomaly is in some ways reminiscent of situations already 
discussed in the literature, where DM is not a thermal 
relic, but rather originates from decays of a heavier state. 
Examples are for instance decays of gravitinos or weakly 
coupled moduli into DM [H^l . Also in these two cases the 
correlation between the annihilation cross section of DM 
and the relic abundance is modified and relaxed from the 
thermal relic relation. In this way also decaying gravitino 
or weakly coupled moduli decays can give large enough 
"boost factors" to explain PAMELA- ATIC anomaly. An 
important difference with 2DM is that neither gravitino 
nor moduli are thermal relics. Rather, their abundance 
(before the decay) is reflective of Planck scale physics. 
The 2DM setup represents the other limit, where X2 in- 
teractions governing its abundance (before the decay) are 
not Planck suppressed and can possibly be probed at fu- 
ture colliders. In this way 2DM model is much closer to 
the simple thermal relic scenario. 

In conclusion, we have presented a simple 2DM model 
that can mimic the large boost factors B needed to ex- 
plain the PAMELA and balloon experiment data. Two 
hierarchies are needed for this proposal to work: i) the ra- 
tio of xi and X2 cross sections need to be large, ~ B, and 
ii) the X2 decay width needs to be much smaller than for 
normal electroweak decays. If the PAMELA data cannot 
be explained using astrophysics, it may be necessary to 
revisit the kind of scenario outlined in this manuscript in 
more detail. 
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